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Abstract 

Ground state energies for plaquette and dimer order in the J1-J2 square- 
lattice spin-half Heisenberg model are compared using series expansion meth- 
ods. We find that these energies are remarkably close to each other at in- 
termediate values of where the model is believed to have a quantum 
disordered ground state. They join smoothly with those obtained from the 
Ising expansions for the 2-sublattice Neel-state at J2/J1 ~ 0.4, suggesting 
a second order transition from a Neel state to a quantum disordered state, 
whereas they cross the energy for the 4-sublattice ordered state at J2/J1 ~ 0.6 
at a large angle, implying a first order transition to the 4-sublattice magnetic 
state. The strongest evidence that the plaquette phase is not realized in this 
model comes from the analysis of the series for the singlet and triplet excita- 
tion spectra, which suggest an instability in the plaquette phase. Thus, our 
study supports the recent work of Kotov et al, which presents a strong picture 
for columnar dimer order in this model. We also discuss the striped nature of 
spin correlations in this phase, with substantial resonance all along columns 
of dimer s. 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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I. INTRODUCTION 



There has been considerable study, over the last decade, of the frustrated spin-| square 
lattice Heisenberg antiferromagnet (the "J1-J2 antiferromagnet"). These studies include 
exact diagonalizations on small systems , spin- wave calculations , series expansions 
|0J^, and a field-theoretic large- iV expansions [^]. 

These studies, and others, have provided a substantial body of evidence that the ground 
state of this system, in the region 0.4 ^ J2I J\ ~ 0.6, has no long-range magnetic order and 
has a gap to spin excitations. For J2I J\ ~ 0.4 the model has conventional antiferromagnetic 
Neel order whereas for J2I J\ ~ 0.6 the system orders in a columnar (vr, 0) phase. Whether 
this "intermediate phase" is a spatially homogeneous spin-liquid, or whether it has some 
type of spontaneously broken symmetry leading to a more subtle type of long-range order, 
has not been conclusively established. 

Zhitomirsky and Ueda |]TU[ have proposed a plaquette resonating valence bond (RVB) 
phase, which breaks translational symmetry along both x and y axes, but preserves the 
symmetry of interchange of the two axes. The horizontal and vertical dimers resonate 
within a plaquette. An early series study [|1T| had investigated the relative stability of 



various spontaneously dimerized states and had concluded that a columnar dimerized phase 
was the most promising candidate for the intermediate region, in agreement with the large- 
N expansions. Zhitomirsky and Ueda claim their plaquette phase has a lower energy 
than this columnar dimer phase, but we find this to be incorrect. 

Further support for the columnar dimer scenario comes from recent work of Kotov et al. 



1^ , who combine an analytic many body theory with extended series and diagonalization 
results to study the nature and stability of the excitations in the intermediate region. It is 
argued that where the Neel phase becomes unstable the system will develop not only a gap 
for triplet excitations but also a gapped low-energy singlet which reflects the spontaneous 
symmetry breaking. This is clearly seen in the calculations. At J2/J1 — 0.38 a second order 
transition occurs, with the energies of Neel phase and dimerized phase joining smoothly, and 
the energy gap and dimerization vanishing. 

It is the aim of this paper to further investigate, using series methods, the competing 
possibilities of columnar dimerization versus plaquette order in the intermediate region of 
the J1-J2 antiferromagnet. It is conceivable that both occur, with a transition from one 
to the other. However, such a transition, reflecting a change of symmetry, is expected to 
be first-order and not well suited to series methods. If both phases are locally stable the 
most direct way to compare them is by comparison of the ground state energies. If one is 
unstable this should show up by the closing of an appropriate gap or by the divergence of 
an appropriate susceptibility. In this paper we calculate the ground state energy and the 
singlet and triplet excitation spectra by series expansions about a disconnected plaquette 
Hamiltonian. We also calculate the susceptibility for the dimer phase to break translational 
symmetry in the direction perpendicular to the dimers. This susceptibility will be large if 
there is substantial resonance in the dimer phase and will diverge if there is an instability 
to the plaquette RVB phase. 

Combining the plaquette expansion results with the dimer expansions of Kotov et al. 



I^ , a very interesting picture emerges for the quantum disordered phase. We find that 



the plaquette phase is unstable and hence is not the ground state for this model. The 
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dimer phase, on the other hand, is stable. However, there is substantial resonance in the 
dimer phase. The spin-spin correlations are not simply those of isolated dimers. Instead, 
the nearest neighbor correlations are nearly identical along the rungs and chains of dimer 
columns. In contrast, the correlations from one dimer column to the next are much weaker. 
The spin-gap phase appears separated from the Neel phase by a second order transition, 
whereas it is separated from the columnar phase by a first order transition. These results 
are in remarkable agreement with the large-N theories The existence of a quantum 
critical point separating an antiferromagnetic phase and a quantum disordered phase with 
striped correlations in a microscopic model makes this critical point a particularly interesting 
one. The role of doping and its implications for high-T^ materials deserves further attention. 



II. SERIES EXPANSIONS AND RESULTS 

We study the Hamiltonian 

H = j,Y.s,-s, + j2 E s^-s,- (1) 

n.n. n.n.n. 

where the first sum runs over the nearest neighbor and the second over the second nearest 
neighbor spin pairs of the square-lattice. We denote the ratio of couplings as y = J2/Ji- 
The linked-cluster expansion method has been previously reviewed in several articles |13|-[15|, 
and will not be repeated here. To carry out the series expansion about the disconnected- 
plaquette state for this system, we take the interactions denoted by the thick solid and 
dashed bonds in Fig. 1 as the unperturbed Hamiltonian, and the rest of the interactions as 
a perturbation. That is, we define the following Hamiltonian 

H = Ho + H, (2) 

where the unperturbed Hamiltonian (Hq) and perturbation (Hi) are 

Hq = Ji E Si ■ Sj + J2 E Sj • Sj 

(ij)eA {ij)eB 

(3) 

Hi = AJi E Si ■ Sj + AJ2 E Sj • Sj 

{ij)ec {ij)eD 

and the summations are over intra-plaquette nearest-neighbor bonds (A), intra-plaquette 
second nearest-neighbor bonds (B), inter-plaquette nearest-neighbor bonds (C), inter- 
plaquette second nearest-neighbor bonds (D), shown in Fig. 1. With this Hamiltonian, 
one can carry out an expansion in powers of A, and at A = 1 one recovers the original 
Hamiltonian in Eq. (p. Thus, although we expand about a particular state, i.e. a plaquette 
state, our results at A = 1 describe the original system without broken symmetries, provided 
no intervening singularity is present. Such perturbation expansions about an unperturbed 
plaquette Hamiltonian have been used previously to study Heisenberg models for CaV409 

It is instructive to consider the states of an isolated plaquette. There are two singlet 
states, one with energy (—2 + y/2)Ji and the other with energy (— 3?//2) Ji. The former is 
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the ground state for y < 1 and corresponds to pair singlets resonating between the vertical 
and horizontal bonds of the plaquette. It is even under a 7r/2 rotation. The latter is the 
ground state for y > 1 and is odd under a 7r/2 rotation. The wavefunctions for these two 
singlet states are 
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where o — o means these two spins form a singlet. There are three triplet states, one with 
energy (—1 + y/2)Ji and a degenerate pair with energy (— ?//2)Ji; like the singlets, these 
have a level crossing at y = 1. Under a 7r/2 rotation the former is odd, while the latter two 
are even and odd, respectively. Finally there is a quintuplet state at (1 + y/2)Ji, which is 
even under a 7r/2 rotation. For y < 1/2 and y > 2 the first excited state of the plaquette 
is a triplet, while for l/2<y<2itis the other singlet. These states and corresponding 
energies are shown in Figure 2. The eigenstates of Hq, the unperturbed Hamiltonian, are 
direct products of these plaquette states. 

To derive the plaquette expansions we identify each plaquette as a 16 state quantum 
object, and these lie at the sites of a square lattice with spacing 2a, where a is the original 
lattice spacing. Interactions between plaquettes connect first and second-neighbor sites on 
this new lattice. The cluster data is thus identical to that used by us previously |^ to derive 
Ising expansions for this model. Because there are 16 states at each cluster site, the vector 
space grows very rapidly with the number of sites and thus limits the maximum attainable 
order for plaquette expansions to considerably less than can be achieved for dimer or Ising 
expansions. 

We have computed the ground state energy Eq to order A^, for fixed values of the coupling 
ratio y. The series are analysed using integrated differential approximants [|l7l, evaluated 
at A = 1 to give the ground state energy of the original Hamiltonian. The estimates with 
error bars representing confidence limits, are shown in Figure 3. For comparison we also 
show previous results obtained from Ising expansions ^ and dimer expansions |]T2|]. We 
find that, in the intermediate region, the ground state energy for both plaquette and dimer 
phase are very close to each other and cannot be used to distinguish between them. The 
dimer expansion yields slightly lower energies near the transition to the Neel phase. We do 
not draw any conclusions from this. 



Zhitomirsky and Ueda |T0[ have claimed that the ground state energy from a second- 
order plaquette expansion is -0.63 (at y = 0.5), much lower that the dimer expansion result 
-0.492. This result appears incorrect. At J2/J1 = \ the ground state energy is given by 



AEo/NJi = -7/4 - 277AV1440 - 0.001357A^ - 0.0210609A^ 

-0.000319586A^ - 0.00580643A^ - 0.001822686A^ + O(A^) 



(6) 
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The second order result (at A = 1) is Eq/NJi = —0.485, rather than —0.63. We note that if 
the second order coefficient were 4 times larger then the resulting energy would be —0.62986. 
We have also derived series, to order A^, for the singlet and triplet excitation energies. 



^sikx,ky), At{kx,ky) using the method of Gelfand fl^, and taking as unperturbed eigen- 
functions the corresponding plaquette states. The low order terms for J2/ Ji = 0.5 are given 
by: 

Asih, ky)/Ji = 1 - 301AV1440 + 137AV86400 + 217A^ cos(A;^) cos(A;j,)/172800 

+ (-5AV16 - 89AV9600)[cos(A;^) + cos{ky)]/2 (7) 

Mkx, ky)lJi = 1 - 3691AV30240 + (-2A/3 + 11AV720)[cos(A;^) + cos(A;j,)]/2 
-A^[cos(2A;^) + cos(2fcj^)]/120 + (A/3 - 5AV96) cos(fc^) cos(fcy) 
-A^[cos(2A;^) cos{ky) + cos{k^) cos(2A;j,)]/90 + 7A2 cos(2A;^.) cos(2A;j^)/360 (8) 

The full series are available on request. We ffist consider the triplet excitations. Figure 
4 shows At{kx,ky) along high symmetry directions in the Brillouin zone for A = 0.5 and 
various coupling ratios y. For A ^ 0.6 the series are well converged and direct summation 
and integrated differential approximants give essentially identical results. We find that the 
minimum gap occurs at (0,0) for J2/J1 ^ 0.55 and moves to (vr, 0) for J2/J1 ^ 0.55. Next 
we seek to locate the critical point Ac where the triplet gap vanishes. This is done using Dlog 
Fade approximants to the gap series at the appropriate {k^., ky). In practice this works well 
when the minimum gap lies at (0, 0). For J2 = we find a critical point at Ac = 0.555(10). 
We can compare this result with recent work of Koga et al. ||18| who obtain Ac = 0.112 
from a modified spin-wave theory and Ac — 0.54 from a 4th order plaquette expansion. The 
critical point Ac increases with increasing y. At y = 0.5, at the approximate centre of the 
intermediate phase, we find Ac — 0.89(7). This result has some uncertainty but, if accurate, 
means that the plaquette phase becomes unstable before the full Hamiltonian (A = 1) is 
reached. The associated critical exponent ly describing the vanishing of the triplet gap is 
about 0.7 for J2/J1 < 0.4, suggesting that the transition lies in the universality class of 
the classical d = 3 Heisenberg model. On the other hand, for J2/J1 ^ 0.4 the exponent 
u is about 0.4. This supports the existence of an intermediate phase lying in a different 
universality class. 

Figure 5 shows the singlet excitation energy As{kx,ky) along high symmetry directions 
in the Brillouin zone for A = 0.5 and the same coupling ratios y as Figure 4. Again the 
series are well converged and direct summation and integrated differential approximants give 
essentially identical results. We find that the minimum gap occurs at (0, 0) for all J2/ Ji- We 
have also noted that for J2/J1 = 0.5, the triplet excitation and the singlet excitation have 
same gap at A = 0, but at A = 0.5, the singlet gap is considerable larger than the triplet gap, 
this means probably that the triplet gap close before the singlet gap at J2/J1 = 0.5. The 
critical point obtained by the Dlog Fade approximant to the singlet gap is also generally 
slightly larger than that obtained from the triplet gap around J2I Ji = 0.5 (see Fig. 6). 

The full phase diagram in the parameter space of J2I J\ and A could be very interesting 
from the point of view of quantum phase transitions, but may not be easy to determine 
by numerical methods. Some possible scenarios are shown in Fig. 7. One possibility is 
that the plaquette phase, for all J2/J11 has an instability to some magnetic phase and the 
dimerized phase exists only very close to A = 1 inside the magnetic phases. A second 
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possibility is that the plaquette-Neel critical line meets the Neel-dimer critical line at some 
multicritical point at a value of J2/ Ji around 0.5, after which there is a first order transition 
between the plaquette and the dimer phases. A third possibility is that the plaquette-Neel, 
Neel-dimer and plaquette-columnar critical lines all meet at some multicritical point. The 
numerically determined phase diagram is particularly uncertain in the interesting region, 
0.5 ^ J2/ Ji ~ 0.6. , where incommensurate correlations could also become important. 

Lastly we have derived expansions for a number of generalized susceptibilities. These 
are defined by adding an appropriate field term 

AH = hY,Q,, (9) 

to the Hamiltonian and computing the susceptibility from 

,, = -l„,„^ (10) 

A divergence of any susceptibility signals an instability of that phase with respect to the 
particular type of order incorporated in x- 

We have computed two different susceptibihties from the plaquette expansion. One is 
the antiferromagnetic (Neel) susceptibility with the operator Qij 

Q,, = i-iy'-'si^ (11) 

The other is the dimerization susceptibility with the operator Qij 

which breaks the symmetry of interchange of x and y axes. We have computed series 
to order for the antiferromagnetic susceptibility and to order for the dimerization 
susceptibility. The series have been analyzed by Dlog Pade approximants. The series for the 
antiferromagnetic susceptibility shows the same critical points (within error bars) as those 
obtained from the triplet gap for J2/J1 ^ 0.4. The series for the dimerization susceptibihty 
is very irregular, and does not yield useful results. For example, for J2/J1 = 0.5, the series 
is: 

Xd = 629/90 + lOlA/300 + 2.0097647A2 - 0.269629A^ + 0.438527A^ + O(A^) (13) 

For completeness, we also compute the susceptibility for the dimer phase to become 
unstable to the plaquette phase from an expansion about isolated columnar dimers [jl2|, by 
adding the following field term: 

AH = hJ2i-iyS,,-S,,+, (14) 

id 

which breaks the translational symmetry in the direction perpendicular to the dimers. The 
series has been computed up to order A^, (note that A here is the parameter of dimerization). 

An analysis of the series shows that this susceptibility becomes very large as A — 1, 
for all J2/J1 and the critical A, where the susceptibility appears to diverge, approaches 
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unity from above as J2/J1 is increased to 0.5. This implies that there are staggered bond 
correlations in the direction perpendicular to the dimers, which extend over a substantial 
range. An interesting question is, in the absence of the plaquette phase as discussed earlier, 
what could these correlations represent? At this stage it is useful to recall another calculation 
by Kotov et al. [|T2|. Within the dimer expansion, they calculated two different dimer order 
parameters. 



and. 



Dx — I < Sjj- ■ Sj+ij > — < Sj+ij ■ Sj-|-2j > I; (15) 



Dy — I < Sjj- ■ Sj+ij > — < Sj.j • Sjj-|_i > I, (16) 



where the elementary dimers connect spins at z, j and i + 1, j. They found that for 0.4 ^ 
J2/ Ji ~ 0.5, Dy is nearly zero, whereas only goes to zero at the critical point. These 
results suggest that the dimer phase consists of strongly correlated two-chain ladders, which 
are then weakly correlated from one ladder to next. This striped nature of spin correlations 
in the dimer phase has not been noted before and is clearly a very interesting result. The 
situation for J2/ Ji ^ 0.5 is again less clear. As discussed before, there are many possibilities 
for the phase diagram in that region and much longer series are needed to throw more light 
on the situation. Perhaps there is an interesting multicritical point in that region of the 
phase diagram. 



III. DISCUSSION 

We have attempted to further elucidate the nature of the intermediate, magnetically 
disordered, phase of the spin-^ J1-J2 Heisenberg antiferromagnet on the square lattice. This 
phase is believed to occur in the range 0.4 < J2/ Ji ~ 0.6. Our approach has been to 
derive perturbation expansions (up to order A^) for the ground state energy, singlet and 
triplet excitation energies, and various susceptibilities, starting from a system of decoupled 
plaquettes (A = 0) and extrapolating to the homogeneous lattice (A = 1). We have also 
derived expansions about an unperturbed state of isolated columnar dimers ("dimer phase"). 
Both of these have been proposed as candidates for the intermediate phase. 

We find that the ground state energy for both plaquette and dimer phases are very 
similar, any difference lying within the error bars. From this result alone we cannot favor 
one phase over the other. 

The analysis of the singlet and triplet excitation spectra suggests an instability in the 
plaquette phase. In particular, in the disconnected plaquette expansions, Dlog Pade analysis 
indicates that the gaps would vanish for A less than unity. The gap appears to close first for 
the triplets and then for the singlets. This is the strongest evidence that the plaquette phase 
is not realized in this model. However, we should mention here that the critical exponents 
associated with the vanishing of the gaps are rather small (< 0.4) and the gap closes not 
too far from A equal to unity. Thus, with a relatively short series, this should be treated 
with some caution. One could ask why the energy series appear to converge well despite 
the instability. However, this is a well known feature of series expansions, that quantities 
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having weak singularities may continue to show reasonable values even if extrapolated past 
the singularity. 

A consistent interpretation of these results is that within the parameter space of our 
non-uniform Hamiltonian, the plaquette phase is first unstable to a magnetic phase, which 
then must give way to the columnar dimer phase. Similar results for the instability of the 
staggered dimer phase were suggested before by Gelfand et al [0. However, the full phase 
diagram in the J2I J\ and A parameter space is difficult to obtain reliably, especially near the 
transition to the columnar phase. There are possibilities of some novel multicritical points, 
which deserve further attention. 

One of our most interesting results is the finding of striped spin correlations in the dimer 
phase. In this phase, the nearest neighbor spin correlations are nearly equal along the 
rungs and along the chains of a two spin column and there are extended bond correlations 
along the chains. However, spin correlations from one column to the next are much weaker. 
In other words, the dimers are strongly resonating along vertical columns. The existence 
of a quantum critical point separating an antiferromagnetic phase with such a quantum 
disordered phase with striped correlations is a very interesting feature of this model which 
deserves further attention in the context of high-Tc materials. 
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FIG. 1. The generalized J1-J2 Heisenberg model with plaquette structure, the couplings Ji, J2, 
xJi, and XJ2 indicated by thick solid bonds, thick dashed bonds, thin solid bonds and thin dashed 
bonds, respectively. 
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FIG. 2. The energies of an isolated plaquette as function of y. The notation presents the 
S'-value and effect of 7r/2 rotation (it) for each state. 
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FIG. 3. The ground-state energy per site Eq/NJi as function of J2/ Ji, obtained from the Ising 
expansions (triangle points connect by dotted Unes), the columnar dimer expansions |jl2| (open 
square points connected by dashed line), and the plaquette expansion (full circle points connected 
by solid line). 




FIG. 4. Plot of the triplet excitation spectrum A(kx,ky) along high-symmetry cuts through 
the Brillouin zone for the system with coupling ratios y = 0, 0.3, 0.5, 0.7 and A = 0.5 [shown in 
the figure from the top to the bottom at {tt,tt), respectively]; the lines are the estimates by direct 
sum to the series, and the points with error bar are the estimates of the integrated differential 
approximants to the series. 




FIG. 5. Plot of the singlet excitation spectrum A(kx,ky) along high-symmetry cuts through 
the Brillouin zone for the system with coupling ratios y = 0, 0.3, 0.5, 0.7 and A = 0.5 [shown in the 
figure from the top to the bottom, respectively]; the lines are the estimates by direct sum to the 
series, and the points with error bar are the estimates of the integrated differential approximants 
to the series. 




FIG. 6. Phase diagram for generalized J1-J2 Heisenberg model with plaquette structure, as 

determined from the plaquette expansions. The full (open) points with error bars and a solid 
(dotted) line to guide the eye indicate the line where the (0, 0) triplet (singlet) gap vanishes. 
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FIG. 7. Some possible topologies for the phase diagram in the A-y plane. Here, N is the Neel 
phase, C is the columnar phase, D is the dimer phase and P is the disordered Plaquette phase with 
no long range order. Note that some of the phase boundaries such as those between Plaquette and 
Dimer and between Dimer and Columnar phases could be first order, whereas the others could be 
second order. M is a possible multicritical point, where several critical lines meet. 
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